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We consider the production of gravitons in an inflationary cosmology by ap- 
proximating each epoch of change in the equation of state as sudden, from which a 
simple analytic graviton mode function has been derived. We use this mode func- 
tion to compute the graviton spectral energy density and the tensor-induced cosmic 
microwave background anisotropy. The results are then compared to the numerical 
calculations which incorporate a smooth radiation-matter phase transition. We find 
that the sudden approximation is a fairly good method. Besides, in determining the 
frequency range and amplitude of the mode function, we introduce a pre-inflationary 
radiation-dominated epoch and use a physically sensible regularization method. 
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I. Introduction 



The spatial flatness and homogeneity of the present Universe strongly suggest that a 
period of de Sitter expansion (inflation) had occurred in the early Universe During 
inflation quantum fluctuations of the inflaton field may give rise to energy density (scalar) 
perturbations 0, which can serve as the seeds for structure formation of the Universe. Also, 
a well-defined spectrum of gravitational waves could be produced from the de Sitter vacuum 
0. Since the waves decouple from matter very early, they remain as a stochastic back- 
ground of gravitational waves at present and provide a potentially important probe of the 
early time. Unfortunately, the detection of these primordial waves by using terrestrial wave 
detectors or the timing of millisecond pulsars (this is only sensitive to short-wavelength 
waves of period less than order years) is still several orders of magnitude below the present 
experimental sensitivity. However, being similar to energy density perturbations, horizon- 
sized gravitational waves can induce distortions of the cosmic microwave background (CMB) 
P^TTI via the Sachs- Wolfe effect ||r^. In fact, large-angular-scale temperature anisotropy of 



the CMB has been recently detected []T3[. So, it is worthwhile to reexamine in detail the 
graviton production from inflation. 

Production of gravitational waves from quantum fluctuations during inflationary period 
have been considered ||^,p!^-|T7| . In transverse traceless (TT) gauge, a classical gravitational 



wave has two independent polarization states. As noticed by Grishchuk |T^, each polar- 
ization state of the wave behaves as a minimally coupled, massless, real scalar field, with 
a normalization factor of vTGttG linking the two cases (we will show this in Section III). 
Thus, the study of graviton production during inflation reduces to considering the quantum 
fluctuations of a scalar field in the de Sitter space-time. The result is the well-known scale- 
invariant spectrum. The specrum can be derived by relating the scalar quantum-mechanical 
two-point function defined in the de Sitter background to a two-point statistical average of 
an ensemble of classical gravitational fields In evaluating the two-point function, the de 
Sitter invariant vacuum has been chosen. The subsequent evolution of the fields is then gov- 
erned by the classical equation of motion. The spectrum can also be derived by a sequence of 
Bogoliubov transformations between creation and annihilation operators defined in various 
stages of the Universe: inflationary, radiation-, and matter-dominated. This consideration 
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reveals that the results are independent of state over a wide range of initial states in the 



inflationary period, thus allowing one to select the de Sitter invariant vacuum W^ . 

In this paper, we will reexamine in detail the graviton production during inflation. We as- 
sume four sequential phases of the Universe: radiation-, inflationary, radiation-, and matter- 
dominated. Each phase change is approximated as a sudden process (i.e., assuming the 
transition from one phase to another phase to be instantaneous). We begin with a scalar 
field in a radiation-dominated period prior to inflation, and work out the dynamics of the 
growing quantum fluctuations of the field during inflation. In evaluating the scalar two- 
point function, we avoid using any formal renormalization schemes to remove divergences, 
instead, we will adopt a physically sensible method introduced by Vilenkin to regularize 
the infinities [|T^. This allows us to determine the frequency range and amplitude of the 
gravitons, and also show explicitly the validity of selecting the de Sitter invariant vacuum 
during the inflationary period. Based on the sudden approximation of the phase changes, 
we will work out the graviton mode function for a wide range of frequency which is relevant 
to astrophysical measurements. This mode function is then applied to compute the spectral 
energy density and the tensor-induced CMB anisotropy. Since in a realistic phase transition 
the phase change is rather smooth than abrupt, we have to examine how good the sudden 
approximation can reproduce the actual situation. In order to do so we will compare the 
results obtained by taking the sudden approximation to the numerical calculations which 
incorporate a smooth phase transition. 

II. Background Metric of the Early Universe 

Let us suppose an initially radiation-dominated universe. The Universe expanded and 
became dominated by the energy density resided in the vacuum which drove the inflation. 
After inflation, the Universe was reheated by the latent heat released from the vacuum and 
resumed radiation-dominated. Later, matter dominated the energy density and the Universe 
is currently in the matter-dominated phase. This supposition makes sense only if inflation did 
occur below the Planck scale (for instance, grand unified theories (GUT) inspired inflationary 
models), since a temperature above the Planck scale may not be defined. At temperatures 
well below the Planck scale but above the GUT scale, radiation most likely dominates the 



energy density of the Universe. For Planck-scale inflation (e.g., chaotic inflation ||20|) or 



eternal inflation (inflationary period extends to negative infinity) , it reduces to a three-step 
calculation in which the de Sitter invariant vacuum in the inflationary epoch is reasonably 
singled out. In the following, we will assume an inflation which is an exponential spatial 
expansion. In fact, our results can be easily extended to other possibiliy like power-law 
inflation. Since the curvature term is sub-dominated at very early times and the Universe 
is essentially spatially flat after inflation, the space-time is well depicted by the A; = 
Robertson- Walker metric, 

ds'' = de - a^{t)d^^ 

{drf - dx^) , (1) 



a 



where airj) and dt] = dt/a(t) are the scale factor and conformal time respectively. Here we 
are using the signature (-1,-1,-1,1). 

Let us consider four phases: pre-inflation radiation-dominated, inflationary, radiation-, 
and matter-dominated. The scale factor for each phase is given by a{rj) — A{ri + 5)^/'^^+^*^ 
according to the equation of state p — qp oi that phase, where A and B are constants. 
For our purposes, we normalize a{rj) = 1 and choose rj — —1/ H at t — when inflation 
begins, where H is the Hubble constant during inflation. Here we will assume sudden phase 
transitions. By requiring that a{r]) and d{r}) (where dot means taking derivative with respect 
to 77) are continuous functions in the cosmological history, we obtain 



a{r]) 



(2) 

-773(77 - 2772)777! , V2<V <Vi 



. -V3{v + Vi- 4772) 74^7! (?7i - 2772) , 771 < 77 
where 7/3,7/2, and 7/1 are respectively the conformal times at which the inflationary, radiation- 
, and matter-dominated eras begin. We have 7/3 = —1/H and 772 = — l/(i/e^^) (where r is 
the duration of inflation). Note that 77 e (2773, -|-oo) and today will be denoted by 770. It can 
be easily shown that 



a(7/i) 

0(772) 

provided that 0(771)70(772), 0(770)70(771) >> 1. 



0(7/0) 
0(771) 
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VI , (3) 



III. Gravitational Wave Equation 



Weak gravitational waves on a Robertson- Walker universe has been investigated by Lif- 



shitz [21|. As we restrict ourselves to a conformally fiat universe, this conformal property 
allows us to rederive the wave equation in a much faster way as follows. 
The theory of pure gravity is given by the Einstein-Hilbert action 



^—jd^x^R, (4) 

TXLr J 
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where G is the gravitational constant. In the weak field approximation, small metric per- 
turbations are ripples on the background metric: 

9i„y = a'^iv){Vi^u + h^^) , hf,^«l, (5) 

where 7]^^, is the Minkowski metric, and Greek indices run from to 3. In this section, we 
will leave a{ri) as a dimensionless and arbitrary scale function. Therefore, without any loss 
of generality, we choose 0(77) very near to unity and write it as 

a{r]) = e"('') ; a{r]) « 1 . (6) 

By expanding Eq. (Q) and keeping terms up to quadratic in h^^, and a, we obtain 

- ]fah''^d^hp^, - 2d^ad^h + 2da(TdpK'^ - Qd^^ad^'a) , (7) 

where all indices are lowered and raised with rj^^, and h = rj^^h^y. In synchronous gauge, 
^00 = hoi = 0, where i runs from 1 to 3. The remaining hij contain a transverse, traceless 
tensor which corresponds to a gravitational wave. Here we are only interested in this tensor 
mode. Henceforth we will work in the TT gauge, i.e., h'^ = dih^^ = and denote the two 
independent polarization states of the wave as +, x \22\. Then, we obtain from Eq. the 
action of graviton as 

1 /• o 1 



Etraviton 



16nG 

We can write a monochromatic wave with a wave vector k as 



J d^x a" -d^hijd^h'^ . (8) 
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hij{x) = h{x; k, A) ey(k; A) , (9) 

where ejj(k; A) is the polarization tensor and A = +, x. The polarization tensor satisfies 

e,,(k;A)e*^(k;A') = 25;,v. (10) 

Hence, for this wave Eq. (j^) becomes 

Won = / d^x a\r]) ^ k, + {d^h{x; k, x))^] , (11) 

which is in fact the action for two real, massless free scalar fields (j){x] k. A) = 
(167rG')^^/^/i(x; k. A) in the background space-time. As a result, each polarizaion state of 
the wave behaves as a real, massless, minimallly coupled scalar field, with a normalization 
factor vTGttG relating the two cases. Writing 

h{x; k. A) = {2n)-hxiv; k) e^^"" + h.c. , (12) 
one finds from Eq. (pUD that the equation of motion for the wave amplitude hx{ri; k) is 

hx + 2-hx + k^hx = , (13) 
a 

where the dot denotes taking differentiation with respect to t]. This is just the Klein-Gordon 
equation for a massless plane wave in the background space-time |T^. In what follows we 
will simply consider a massless real scalar field 0(x). 

IV. Solutions of Wave Equation 

We consider the plane-wave mode expansion of a massless real scalar field in the back- 
ground metric (p, 

(pix) = (27r)-t J d^k [a^Mv) e"^"" + h.c 

[ak,aj,,] =5(k-k') , (14) 

where the mode function ipkiv): which defines the vacuum state, satisfies Eq. (|1^). When 
the scale factor a{r]) takes the form given by Eq. (0), we find: 
(i) When r] < rj^, 
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(2), 

2 



(15) 



where ^ = 2/H + r], H^^ and HY' are the Hankel functions. The coefficients ai and 0:2 are 



r(2) 



functions of k, subject to the normahzation condition: 



I |2 I |2 1 



(16) 



(ii) When rj^ < rj < 772, 



/3iM^)(A:r7)+/32M'^(A;r7) 

2 



(2), 

2 



(17) 



(iii) When ri2 < r] < rji, 



Mr^)=a-\r])-{n\^\)-^ 



2 



(2)^ 

1 

2 



(18) 



where ^ = Tj — 2ri2. 
(iv) When 771 < 77, 



^,(7^) = a-i(r^)^(7r|e|)^ 



2 



(2), 

3 

2 



(19) 



where ^ = + ^1 — 4:^2- Similar to the case (i), the Greek coefficients in cases (ii)-(iv) are 
k dependent and each pair is subject to the same normahzation condition (|16|). We hst 
below the Hankel functions and their properties which will be useful for determining the 
coefficients: 



Hl^\-z) = e-''''Hl^\z) - 26-"'"' J^{z) 
Hi^\-z) = e-'''''Hl^\z) + 2e^-V,(^) , 

V TTZ 




2 

i\l — e 



-(l + -)e^ 

TiZ Z 



, i/r(.) = -J-(l--)e-, 

2 \ HZ z 



(20) 



where Jy and are Bessel functions, and are spherical Bessel functions, z > Q 
is real, and v is complex. 

To determine the coefficients, we require that ipk and its time derivative ipk are continuous 
at Tj equal to 773, 772, and r/i. Thus, referring to Eqs. ( p!5| - |T9|) , we have six equations with eight 



unknowns. However, for rj < rjs, we choose the conformal vacuum state, i.e., ai = and 
^2 = 1 in Eq. (0), which corresponds to the positive-frequency mode functions, 'ipkiv) = 
i a^^{7]) e~*'^^/v^2A;, where ^ = 2/11 + 7]. Note that these mode functions are conformal 
transforms of the positive-frequency mode functions defined in Minkowski space-time. As a 
result, the remaining six coefficients can be uniquely determined. In terms of 

Z3 = -kr]3 = T7 ' ^2 = -kr]2 = t?^^ > = HVi - ^V2) = t? o/ x > (21) 



we find that 



2zl ' 



i 1 



/32 = I 1 - - - ^ ^ 



I 1 \ . ^ 

2 / ^ 2 ' 



,, = _,e-,,^-,e-.-„(l + 5i-__L). (22) 

These results are obtained for any arbitrary k. In the next section, by examining the 
dynamics of the growing quantum fluctuations of the scalar field during inflationary period, 
we will show that only plane-wave modes of k within certain range will be excited and 
generated from the de Sitter vacuum. 

V. Scalar Field Fluctuations 

The general formalism for calculating the quantum fluctuation (0^) in the de Sitter space- 
time has been developed by Bunch and Davies p3|. Using a point-splitting regularization 



scheme they calculate a two-point function for in Eq. ([T^), 



= /J ^'^■^'''"'HMWkW^k , (23) 



oo 



where ipkiv)^ which defines the vacuum state, is given by Eq. {^7^. Furthermore, they 
specialize to the de Sitter invariant vacuum such that /5i = and (32 = 1- With this 

8 



choice, Eq. (^3]) is singular for x" — ^ x' . This singularity can be removed by an appropriate 
renormalization scheme. In the present paper, we will treat (0^) based on an alternative 



procedure introduced by Vilenkin |T9 



In an inflationary universe one has to carefully consider the physically relevant quantity 
to be calculated. As has been discussed by Vilenkin, in the inflationary universe the modes 
of interest are only those with wavelength greater than the horizon, as these are those 
responsible for the growth of (0^). These are the modes with k\!]\ << 1, so that the 
integration over k should be cut off at A; ~ He^'^. When k < H, the magnitude of the 
fluctuation depends on the initial conditions of the universe, but with a power-law expansion 
rate of the universe prior to inflation there is no reason to believe that (0^) exceeds the 



anomalous de Sitter fluctuations |T9[. Hence we can write (0 ) at time t as 

(0^) = (0^)o + — / dke\M^)\\ (24) 



where (0^)o is the initial value and t < r. Although Eq. (24) is not without arbitrariness, 
Vilenkin's regularization tecnique makes physical sense as well as being useful from a cal- 
culational point of view. In fact, this method can easily reproduce previous results. For 
example, for a massless scalar field, the growth of quantum fluctuations during inflation in 
this formalism, by using Eqs. (P^D, (^D])-(^), and (P^, is calculated as 

(0') = ^\ri? dk k'\f3,H^P{k7^)+(32Hf\kr^)\' 

—^Ht as t » H-^ . (25) 

This is the well-known result of the linear growth of fluctuations P, p^p4| . In essence, short- 
wavelength quantum fluctuations in the scalar field ipk of wavenumbers, H < k < He^"^, 
will get red-shifted out of the horizon during the inflationary period, after which they freeze 
in, remaining with constant amplitude (see below). The superposition of these frozen modes 
therefore constitutes the coherent scalar field. When the mode re-enters the horizon much 
later during the radiation- or matter-dominated era it appears as a long-wavelength, classical 
wave. 

VI. Graviton Mode Function 
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We write the graviton field in terms of tlie plane-wave modes given in Eqs. (^) and ([l2D 

as 



hij{x) = (27r)-i J2 I [«A(k) hx{Ti] k) e^^'' ey(k; A) + li.c] , 

[aA(k),4,(k')]=5(k-k')(5AA'. (26) 

As explained above, the gravitational wave is related to the scalar field by h\{ri\ k) = 
{IGTiGy^'^ipki'T]). In the background metric (0), ipk{f]) are given by Eqs. (|T5|)-(|T9D with 
the coefficients given by Eqs. (^T])-(^2]). Similar to the scalar case, only quantum fluctua- 
tions in the gravitational field hx{ri; k) of short wavelengths, H < k < He^'^, will get pushed 
outside the horizon during the inflationary period. Therefore, from Eq. (0), z^, Z2, and zi 
are such that 

l<Z3<a{r]2), —-<Z2<1, -4r^<zi<4^- 27 
aim) a (^2) a{ri2) 

In an inflationary cosmology, the exponential expansion factor 0(772) is at least 10^^ in order 
to circumvent the cosmological problems Let us assume aim) = 10"^^. If the reheating 
temperature after inflation is about the GUT scale of order lO^^GeV, then the ratio of 
the scale factors at the beginning of matter-dominated era and at the end of inflation, 
a{rji) / a^m) 1 would be about 10^^. In addition, the present scale factor airjo) ~ 10^a(?7i) p5| . 
For definiteness, we fix 

a(,2) = 10-, 4^^ = 10-, = 10^ (28) 

Hence, we find from Eq. (|27|) that 

1 < ^3 < 10^^ , 10"^^ <Z2<1, 10"^ <zi< 10^^ . (29) 

It is useful to compare the size of a wave to the horizon size at time 77. For the background 
metric (|1]), after inflation, the horizon size / grows as 

I = a(r]) / (i?7 ~ a{r])r] , for 77 >> 772 . (30) 

If we consider a wave with physical wavelength Aphys which is just entering the horizon at 
time f], then its wavenumber k must satisfy 
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k 



2tt a{r]) 2tt a{r]) 2tt 



(31) 



Aphys I V 

This condition separates two limiting cases: a wave with kt] << 27t is well outside the 
horizon whereas a wave with kr] » 2tt is well within the horizon. 

Of particular interest are short-wavelength gravitational waves of period which is within 
the measurable range of terrestrial or astrophysical wave detector, as well as waves of long 
wavelengths comparable to the present horizon size. We can express the wavenumbers k of 
these waves in terms of the present time r/o by using Eq. (|31|) as 

2mi 



k 



(32) 



where n = 1 corresponds to a wave just entering the present horizon, and the wave with 
n > 10^^ (this wave re-enters the horizon during the radiation-dominated epoch at time 
7] < lO^^'^r^o) has present period less than about a year. Substituting Eq. ( P^ in Eq. 
and then using Eqs. (^) and (^), we find for such waves that 



Z3 ~ lO'^nvr , ^2 — 10 ^^nvr 



lO^^nvr . 



(33) 



Also, Eq. (p9D implies 10 < n < 10 . Therefore, in either case, Z3 >> 1 and Z2 « 1. 



This leads to approximation of expressions (|22|) as 



7i - 72 
61 ~ 7i 

52 ~ 7i 



t e 



-I 

1 



+ i e 



— z e' 



3i2i 



8^2 
1 " 

8if 



(34) 



where z^, z^ and Z\ are given in Eq. (^). By using this and Eqs. (|T^)-(|T9|), we can construct 
the wave form of the gravitational wave for any particular value of n. For example, at time 
rj » r]i in the matter-dominated epoch, the wave amplitude is found to be given by 



hxiv, k)k 



3 

2 ~ 



(327rG)^if^e-2*^^ 



kr] 



% e 



2^1 ^zl 



}ii\kr]) 



I 1 + 



% 1 

2^ " Sif 
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I e 



/if^(A;r/) 



\zi 



(35) 



where we have used ^ = r] + rii—4r]2 ^ f] in Eq. (0) since rj >> rji >> 772- In next section we 
will use this expression to calculate the effects of primordial gravitational waves on the CMB. 
We note that Eq. ( ^5]) is at variance with previous results as found in Refs. However, 
we will show shortly that those results can be in fact reproduced by taking different limits 
of Eq. (^). Also, similar expressions can be found in Refs. P^ , pr5|pB| , p7| . In the radiation- 
dominated epoch, for 7] >> ?72 5 we find that 

hxiv, k)kl ^ -{8nG)^He-^'''jo{k7]) , (36) 



which agrees with the result given in Refs. [p!4tp|. 

As mentioned above, a gravitational wave well outside the horizon will remain with 
constant amplitude until much later it re-enters the horizon as a propagating classical wave. 
To illustrate this, we consider a gravitational wave with wavelength larger than the present 
horizon size {n < 1). For this wave, zi < lO^^vr << 1. Thus, further approximation of 
Eq. ( ^4|) can be made: 

3 

5i~(52~7i— , (37) 

where 71 is given in Eq. (0). By using Eqs. ([19|) and (0), we obtain the amplitude squared 
of the wave in the matter-dominated epoch as 

\h,{^; k)\' ^ 16nGa-\v)^A\5,\'^f,m , (38) 

4 TT 



where ^ = tj + rji — 4?72. For rj » rji, this can be simplified as 



\hxiri;k)\^k^ ^ 8nGH 



3ji{kri) 



1 2 



krj 



(39) 



where 

3ji(A;r/) 



^ 1 as n << 1 . (40) 



krj 

The expression (^9]) which is the long-wavelength limit of Eqs. (^9D or (^5]) coincides with 
the well-known result of the scale-invariant spectrum of gravitational waves generated by 



inflation |@-|S,14]. After a wave has entered into the horizon, its wave amplitude decreases 
with time. To show this, we consider a wave with n = 10^°, which corresponds to zi = 
lO^TT >> 1. Then, for rj >> rji, we find from Eq. (|35|) that 
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\hx{r];k)\^k^ ~ 12SnGH\kr] 



(41) 



krj 

which agrees with the result given in Ref. |7[] except differing by a factor of 16. Note also 



that the /c-dependence in Eq. (^) differs from that in Eq. (^), namely, the former has 
less power suppression in k. This behavior has also been found by using transfer function 
methods 



VII. Sudden Transition Approximation 

The derivation of the graviton mode functions in Sec. VI is based on the sudden approx- 
imation of the phase changes. One should justify this approximation before he can trust the 
result in Eq. (0). As seen in Eq. (Q), we have approximated each phase transition as being 
instantaneous. In reality, changing from one phase to another phase must take some finite 
time. The phase transition from the pre-infiationary radiation-dominated era to the infla- 
tionary era and the reheating at the end of inflation should involve very short time scales, 
hence only affecting the extremely high frequency part of the graviton spectrum. Discussion 
of the influence of the reheating on the high frequency cutoff of the spectrum can be found 
in Ref. p9[. Here we will mainly concern the affect of the radiation- matter phase transition 
on the comparatively much lower frequency region of the spectrum. We will see that this 
phase transition would affect the spectrum in the frequency range that is relevant for CMB 
anisotropy calculation. 

Many attempts have been made in the past in finding an accurate solution for the graviton 
mode function across the radiation-matter transition. In a two-component (radiation plus 
dust) cosmology which has a smooth radiation-matter phase transition, the equation of 



motion (|T^) can be solved exactly in terms of complicated spheroidal wave functions |T6| , |30 
In addition, the equation can be solved approximately by the time-independent or time- 
dependent |^l| transfer function method, or by the WKB approximation |3^. Of course, the 



equation can be easily solved via numerical method. However, sometimes it is much more 
convenient to have a simple yet reasonably accurate graviton mode function. Undoubtedly, 
the sudden approximation is the most economical way to obtain such a simple analytic 
solution. 
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To justify the accuracy of the sudden approximation method, we compare the period 
of a wave, 27r//c, with the characteristic duration of the radiation-matter phase transiton, 
Ai], at time rji. Unfortunately, the phase transition, being rather gradual than sudden, does 
not bear a characteristic timing. However, Arj should be comparable to rji. As such, the 
condition that the phase transition is sudden as seen by the wave is 

2vr , 

— >> Ar] ~ r^i , i.e. , Zi << 2it . (42) 
k 

Therefore, Eq. (|35| ) is valid for zi << 2tt, i.e., for waves entering the horizon during the 
matter-dominated era. However, we should look into this matter in more detail. Let us 
recast Eq. (|13]) in terms of the function g\ = a{r])h\ in the form 

9x+{k'-^^g, = 0, (43) 

where 

^ = ("' """'^ (44) 

for a sudden radiation-matter phase transition. In the radiation-dominated era, the solution 
to Eq. ( ^3|) which matches the wave function in the inflationary phase is given by sm{kvi). 

3 

Hence, hxk^ oc sm.{k'r]) / [ka) (where a oc 77), which is in fact the result in Eq. (pq). In the 
matter-dominated era, the modes whose /c's satisfy kr] > zi >> \f2 will not be affected 
by the term aja. Then, the solution to Eq. (|43|) for these modes is again given by •simikr]). 
Hence, the matter-dominated mode function hx is given by the same solution in the radiation- 
dominated era with a ccrf instead. This actually reproduces the result in Eq. (^T]). Although 
the phase transition as seen by the short-wavelength modes is gradual, they are not affected 
by the sudden phase change at all and only scale as 0(77). This is actually an illustration 
of the adiabatic theorem. Therefore, the sudden approximation can give good approximate 
graviton mode function across the radiation- matter phase transition for modes with Zi « 
2tx 01 Zi » \/2. In the next section, we will use Eq. (|35|) to calculate the present spectral 
energy density of primordial gravitational waves and the induced CMB anisotropy. These 
will be compared to the results obtained by considering an equation of state which smoothly 
interpolates between the radiation- and matter-dominated eras. We will find that the sudden 
approximation can actually give fairly good approximation of these quantities even in the 
regime ~ 1. 
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VIII. Spectral Energy Density and CMB Anisotropy 



To incorporate a smootli radiation-matter phase transition, we consider a two- 
compontent universe containing radiation and dust. For convenience, we define the new 
variables 



k' 



krji 



V2-I' 



(45) 



It can be easily shown that 



(46) 



Note that the radiation-matter equality time ri[ ^ 0.41. For alrj^) / a{ri[) = 6000, the present 
time ?7o — 76.5. (Note that ?7q ~ 63.5 if we have used Eq. (|^).) Then, the gravitational 
wave equation ([1^) can be exactly solved in terms of complicated spheroidal wave functions 
16| , |30[| . Here we simply solve it numerically, using the following initial conditions for hy. 

dhx 



hxk' 5 



-^"'{8nG)-^H and 



dr]' 







as 



V^O, 



(47) 



which is consistent with Eq. (|5B|). 

From Eqs. (0) and (^61), the spectral energy density of gravitational waves can be written 

as 



Pg 



A=+,x "''^ A=+,x 



1 




k"\hx\' + 


dhx 


2" 


32TTGa^{r]') 




dr]' 





(48) 



It is useful to define a dimensionless parameter which is the spectral energy density divided 
by the closure density of the Universe, 



Pc 



StcG 



(49) 



where the Hubble parameter H{ri') = [da/ dr]') / a?' . As the graviton production rate during 
inflation for each polarization state is statistically equal, it is expected that h^ = h^ = 
_^-2tz3 jji^/-3/2 where h is real and normalized to unity as i]' — > 0. Hence, 



^ 97r ' 



n' + l, 



(50) 
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where v = SGif^/Svr = Vo/mp; is the inflation parameter (Vq and mpi are the de Sitter 
vacuum energy and Planck mass respectively). To estimate the amplitude of this quantity, 
one typically makes use of the following physical argument. Before a gravitational wave 
enters the horizon, its amplitude is constant. After it enters the horizon, on the other hand, 
it behaves effectively as radiation and will scale with a as such. This gives a fi^ ^ 10~^^ 
for V ~ 10~^ Since we have already derived the graviton mode function, we can 

directly calculate fig at any time. Fig. 1 shows the present time Qg in units of v verus 
the wavenumber in terms of n. Note again that n = 1 corresponds to present horizon- 
sized waves. The curves denoted respectively by the symbols hg, and are obtained 
by using the numerical solution of Eq. (0) with the initial conditions (0), the sudden 
approximation mode function (|35|), and the matter-dominated mode function (pQ]). We see 
that in the horizon- and super horizon-sized wavelength region, the three curves are equal. 
While hm severely underestimates the spectral energy density at large n's, hg is a fairly good 
approximation for all n's. The worst underestimation by hg is by a factor of 2.5 at n ^ 35. 
At large n's, it approaches to a constant plateau and then falls off at the high frequency 



cutoff P5| , p9| , being consistent with the numerical result. 

In the last section, we have argued that the sudden approximation should be valid for 
any n which is much smaller or larger than about 100 (corresponding to 2:1 ~ vr). This is 
verified in Fig. 1. It has been suggested that 61 and 62 in Eq. ( p4D for zi > 2n (or n > 200) 
should be modified to exponential froms (see Eqs. (4.22) and (4.23) of Ref. ||2^) due to the 
adiabatic theorem. However, we find that the spectral energy density obtained by using the 
sudden approximation (|35|) differs by less than 34% for n > 200 from the numerical result. 

We now turn to calculate the CMB anisotropy induced by the primordial gravitational 
wave background. The temperature anisotropy is induced via the Sachs- Wolfe effect, 

5T I r ■ d 

^(e) = -- I rfAeV— /i,,(V,f), (51) 



where e is the propagation direction of the photon, hij is a field operator given by Eq. 
and the lower (upper) limit of integration represents the point of emission (reception) of the 
photon. In CMB measurements, the measured temperature anisotropy is usually expanded 
in terms of spherical harmonics. 
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ST 

— (e) = ^ aimYim (e) . (52) 

l,m 

Since 6T/T arises from quantum fluctuations generated during inflation, it is a Gaussian 
random fleld. This implies that independent Gaussian random variables satisfying 

{O'lmO'Vm') = 2/ _^ ^^ll'^mm', (53) 

where Ci is the anisotropy power spectrum 0], from which we can construct the two-point 
temperature correlation function, 

(^(ei)^(e2)) = ^ECiPiiei ■ e,), (54) 
T T An I 

where Pi is a Legendre polynomial. The formula for the power spectrum is given by |^ 



1 2 



(55) 



where ji is a spherical Bessel function. We will use the decoupling time r]'^^^ ~ 1.54 corre- 
sponding to a{Vo) / (^iVdec) = 1100. In Fig. 2, we plot C; for / = 2 — 250 by using the three 
different mode functions for h as in the calculation of the spectral energy density. Since the 
time integration in Eq. (^) is getting close to the radiation-matter equality time, we should 
keep C, = 1] + rji in Eqs. (^) and (|39D (i.e., replacing t] in these equations hj r] + rji). We 
see from Fig. 2 that the three curves almost match for / < 20. The matter-dominated mode 
function hm severely underestimates the Ci for large /. Using the sudden approximation 
mode function hs, on the other hand, we obtain much better accuracy; the Ci we flnd differs 
by less than 2% for / < 30, less than 4% for / < 130, less than 10% for / < 150, and less 
than 25% for / < 250 from the numerical result. 

For a flxed I, the main contribution to the integral ( p5D for Ci comes from the mode of 
wavenumber k' c:^ l/rjQ (or / 2nn) at the horizon crossing time rj'^ ^ n/k' PJ3^. In virtue 
of this, the mode functions hg and h^, while underestimating the spectral energy density at 
large n, underestimate the anisotropy power spectrum at large /. It is interesting to note 
that in Fig. 2 all three curves have minima at / ~ 170. While both hn and hs have maxima 
at / ~ 220, hm has a maximum at a lower / ~ 210. It has been pointed out that it is the 
temporal phase of the gravitational wave (the phase of h) which determines the overall shape 
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of the power spectrum at large / . Specifically, the location of a maximum or a minimum 



is controlled by the phase of the wave. We thus plot hg and each of wavenumber 
k' ~ 2.2 corresponding to / ^ 170 at near the horizon crossing time r]'^ ~ 1.4 in Fig. 3. We 
find that both hs and hm-, while underestimating the amplitude, reproduce accurately the 
phase of the wave, in accordance with the behavior of the curves at / ~ 170 in Fig. 2. In 
Fig. 4, we plot again the three mode functions with wavenumber k' ~ 2.9 corresponding to 
/ ^ 220 at near t]'^ ~ 1.1. We find that hg reproduces rather good the phase of the wave and 
the location of the maximum at / ~ 220. But, at this /, not only hm badly underestimates 
both the amplitude of the wave and the power spectrum, but also its Ci curve rises up too 
quickly from the minimum to the maximum. 

The tensor-induced CMB anisotropy power spectrum has been calculated by using nu- 
merical codes to evolve the photon distribution function using the general relativistic Boltz- 
mann equation for radiative transfer [|ll|,Q, and by semi-analytic approaches using the 
Sachs- Wolfe integral (0). To evaluate the integral, different graviton mode functions of 
varying accuracy have been used [^, ^ , ^0| -|32| . We find that the hn curve in Fig. 2 is very 



close to the result from Boltzmann codes |Q as well as the result from spheroidal wave 



functions Note that we have used rj'^^^ ^ 1.54 in Eq. ( |55D for all three cases. But, in the 
sudden approximation, Eq. (||) implies that r/^^^ should have been equal to 1.15. However, in 
order to make comparison between the three curves, we insist in using a common rj'^^^ ~ 1.54. 
If we have used rj'^^^ ~ 1.15 to evaluate the integral (^5]), the hg curve should have been simi- 



lar to the results in [|27| , p0| . There they have also used a similar sudden approximation mode 
function to calculate the power spectrum. Ref. has proposed an approximate graviton 
mode function which in fact has an accuracy comparable to this work, but has to involve 
more complicated time-dependent transfer functions. 

VIII. Discussions and Conclusions 

In this paper, we have attempted to derive a simple yet reasonably accurate analytic 
mode function for the gravitational waves generated in inflationary cosmology. The wave 
amplitude in Eq. ( pSf ) is our main result. To derive this result, firstly, we have assumed 
that the Universe underwent four phases in sequence: pre-infiation radiation-dominated, 
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inflationary, radiation-, and matter-dominated. This is quite natural as long as one con- 
siders GUT-scale inflation. To the first order approximation, we have treated each phase 
transition as being instantaneous. Secondly, we follow Vilenkin's physically sensible method 
to regularize the infinities in the evaluation of the scalar two-point function. In essence, the 
regularization scheme simply introduces an infrared cutoff and an ultraviolet cutoff. As a 
result, only those quantum fluctuations in the gravitational field with wavenumbers lying 
between the two cutoffs get pushed outside the horizon during the inflationary period, after 
which they freeze out and become collective modes. Our four-phase calculation explicitly 
shows that (see Eq. (p2D ) it is too naive to select the de Sitter vacuum = 0,^2 = 1) 
during the inflationary phase. Nevertheless, it could be a fairly good approximation since 
Z3 > 1 (see Eq. (pQ])). Moreover, for gravitational waves which are relevant to CMB observa- 
tions, pulsar timing measurements or terrestrial wave detectors, it is perfectly fine to choose 
the de Sitter vacuum (see Eq. (0)). Besides, we can reproduce previous results by taking 
certain limits of Eq. (|53|). For gravitational waves with large wavelength comparable to 
the present horizon size, we have obtained Eq. (|3^) which is the well-known scale-invariant 
spectrum. For shorter waves, the wave form is different from Eq. (|39|), and only by choosing 
some special value of the wavenumber can one obtain the result given in Eq. (P|). 



We have used the sudden approximation mode function (|35D to calculate the present 
graviton spectral energy density and the induced CMB anisotropy. The results are compared 
to the results obtained by using the numerical solution for the graviton mode function in 
a two-component universe which incorporates a smooth radiation-matter phase transition. 
We have found that the approximation of the phase change as a sudden process is a fairly 
good method. Although the mode function ( P^D mildly underestimates the spectral energy 
density and the anisotropy power spectrum, it reproduces accurately the overall shape of 
the power spectrum. This simple analytic mode function should be quite useful in studying 
tensor perturbations from inflationary models. 
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Captions 



Figure 1. Present graviton spectral energy density in units of v. The symbols hn, hg, 
and hm denote respectively the curves obtained by using the numerical mode function, the 
sudden approximation mode function, and the matter-dominated mode function. Note that 
n = 1 corresponds to a wavelength of the present horizon size. 

Figure 2. Anisotropy power spectra by using hn, hg, and respectively. 

Figure 3. Graviton mode functions each of wavenumber k' ~ 2.2 corresponding to / ~ 170 
at near the horizon crossing time r]'^ ~ 1.4. The curves in descending amplitude are hs, 
and hm respectively. 

Figure 4. Graviton mode functions each of wavenumber /c' ~ 2.9 corresponding to I ~ 220 
at near the horizon crossing time 77^ ~ 1.1. The curves in descending amplitude are hn, hs, 
and hm respectively. 
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